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Studies of black hole formation from gravitational collapse have revealed interesting non-linear 
phenomena at the threshold of black hole formation. In particular, in 1993 Choptuik studied the 
collapse of a massless scalar field with spherical symmetry and found some behaviour, which is quite 
similar to the critical phenomena well-known in Statistical Mechanics and Quantum Field Theory. 
Universality and echoing of the critical solution and power-law scaling of the black hole masses have 
given rise to the name Critical Phenomena in Gravitational Collapse. Choptuik's results were soon 
confirmed both numerically and semi-analytically, and have extended to various other matter fields. 

In this paper, we shall give a brief introduction to this fascinating and relatively new area, and 
provides an updated publication list. An analytical "toy" model of critical collapse is presented, 
and some current investigations are given. 

PACS numbers: 04.70.Bw, 04.25.Dni, 04.40-b, 04.50.-Hh, 98.80. Cq 



I. INTRODUCTION 

Gravitational collapse of a realistic body has been one 
of the most important and thorny subjects in General 
Relativity (GR) since the very early times of GR [Q. 
The collapse generally has four kinds of possible final 
states. The first is simply the halt of the process in a 
self-sustained object, such as, stars. The second is the 
dispersion of the collapsing object and finally leaves be- 
hind a flat spacetime. The third is the formation of black 
holes with outgoing gravitational and matter radiation, 
while the fourth is the formation of naked singularities. 
For the last case, however, the cosmic censorship hypoth- 
esis declares that these naked singularities do not oc- 
cur in Nature. 

Due to the mathematical complexity of the Einstein 
field equations, we are frequently forced to impose some 
symmetries on the concerned system in order to make the 
problem tractable. Spacetimes with spherical symmetry 
are one of the cases. In particular, gravitational collapse 
of a minimally coupled massless scalar field in such space- 
times was studied both analytically ||^ and numerically 
Q, and some fundamental theorems were established. 
Quite recently this problem has further attracted atten- 
tion, due to Choptuik's discovery of critical phenomena 
that were hitherto unknown As a matter of fact, it 
is so attractive that Critical Phenomena in Gravitational 
Collapse has already been a very established sub-area in 
GR, and several comprehensive review articles already 
exist 

In this paper, a summary of an invited talk given at the 
XXI Brazilian National Meeting on Particles and Fields, 
we shall first briefly review the subject and give an up- 
dated list of publication in this area. This will be done 
in Sec. II, while in Sec. Ill, we shall present an analytic 
"toy" model of a collapsing massless scalar field. The 



word "toy" model here means that the model doesn't 
really represent critical collapse, since the perturbations 
of the corresponding "critical" solution have more than 
one unstable mode. However, it does have all the main 
features of critical collapse. Since so far, no any critical 
solution is known explicitly in a close form, this toy model 
still serves as a good illustration to critical phenomena 
in gravitational collapse. The paper is closed by Sec. IV, 
in which some current investigations in this fascinating 
area are given. 



II. CRITICAL PHENOMENA IN 
GRAVITATIONAL COLLAPSE 

Starting with spherical spacetimes, 

ds^ = -a^{t, r)dt^ + a^{t, r)dr^ + r^dfl'^, 



(2.1) 



where dfl = dO^ + siii^ 9 dip"^ , and {x^} = {t,r,d,ip} are 
the usual spherical coordinates, Choptuik Q investigated 
gravitational collapse of a massless scalar field, (j>, which 
satisfies the Einstein-scalar field equations. 



□ = 0, 



(2.2) 



where R^^i, denotes the Ricci tensor, K EE [87rG/c'*] is the 
gravitational couphng constant, ( ),^ — d{ )/dx'^, □ = 
g^^S/aiS/ f3, and Va denotes the covariant derivative. 
Once an initial smooth configuration of the massless 
scalar field is given, these equations uniquely determine 
the later evolution of the spacetime and the scalar field 
Q. Let the initial distribution of the massless scalar 
field be parameterized smoothly by a parameter p that 
characterizes the strength of the initial conditions, such 
that the collapse of the scalar field with the initial data 
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p > p* forms a black hole, while the one with p < p* does 
not. A simple example is the gaussian distribution of the 
massless scalar field 



exp 



r - To 



(2.3) 



where to denotes the initial time of the collapse, and 
(j)o, To, S, and q are constants [See Fig. 1]. 
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FIG. 1. The initial configuration of the massless scalar field 
a,t t — to given by Eq.(2.3) in the text. It actually represents 
a collapsing spherical shell, made of the massless scalar field, 
with its thickness S and centralizing at the radius r = ro. (j)o 
represents the amplitude of the wave packet. 



In this case, Choptuik found that the parameter p can 
be any of the four parameters. 



P = {(f>Q, ro, q, S}, 



(2.4) 



that is, fixing any three of the four parameters, for ex- 
ample, ro, 5 and q, and leaving only 0o change, we shall 
obtain a family of initial data, S'[(/>o]- For this family of 
initial data Choptuik found that there exists a critical 
value (/)o such that when (/)o > </)g the collapse always 
forms black holes, and when (po < (t>o the massless scalar 
field first collapses, then disperses to spacelike infinity, 
and finally leaves a flat spacetime behind without forming 
any kind of spacetime singularities. When (I)q w ^J+A^Oj 
where A(f)o is very small, after certain time all the col- 
lapses are quite similar and approach to the critical one 
{(j>o = <Po)- But, at the very end, the collapse will sud- 
denly runs away from the critical one, by either forming 
black holes or dispersing to infinity, depending on the 
signs ofA^o- Choptuik found that for the configuration 
of Eq. ( |2.3| ) there are four different families of initial data, 
S'[(/)o], S[ro], S[S] and S[q], which all exhibit the above 
behaviour. 

In addition to these four, Choptuik also studied many 
others and found that for all the families that behave 
as above, the so-called generic smooth families of initial 
data, the critical solutions are identical, or in another 



word, universal. Moreover, the critical solution is also 
periodic, that is, 



where A* 



{a* , a* and 



In 



eo(r), 



(2.5) 



(2.6) 



with ro being a dimensionful constant, and ^o(''") a pe- 
riodic otherwise arbitrary function with period A. The 
constant A is a dimensionless constant, which was nu- 
merically determined as A « 3.447. 

Yet, near the critical solution but with p > p* , the 
mass of black holes takes the scaling form 



BH 



K{p~p*r, 



(2.7) 



where K is a family-dependent constant, but 7 is another 
dimensionless universal constant, which was numerically 
determined as 7 « 0.37. 

Universality and echoing of the critical solution and 
power-law scaling of the black hole masses have given 
rise to the name Critical Phenomena in Gravitational 
Collapse. 

Choptuik's results were soon confirmed by several inde- 
pendent studies both numerical jl^ and semi-analytical 
|l5|, and have been extended to other matter fields, such 



as. 



• Axisymmetric gravitational waves [[l6| ; 

• Perfect fluids with the equation of state p — kp, 
where p denotes the pressure of the fluid and p the 
energy density and fc is a constant; 

• Quantum black hole formation in 2-dimensional 
spacetimes 

• Non-linear cr-models in two dimensional target 
space ||20| ; 

• Massless scalar field in Brans-Dicke theory 

• SU{2) Yang-Mills field ||; 

• Einstein-Maxwell-scalar fields |Q ; 

• Massive scalar field p^ ; 

• Gravitationally collapsing primordial density fluc- 
tuations in the radiation dominated phase of the 
early Universe ||25[| ; 

• SU{2) Skyrme fleld |§; 

• The collapse of coUisionless matter of the Einstein- 
Vlasov equations p7| , p8[ ; 

• Topological domain walls interacting with black 
holes \M\ 



2 



The gravitational collapse of massless scalar field 
in higher dimensional spacetimes 30 1; 



• Non-linear cr-models in three dimensional target 
space ||3l| ; 

• Gravitational collapse in Tensor-Multi-Scalar and 
Non-linear Gravity Theories 

• Boson stars ^] ; 

• Massless scalar field coupled with the cosmological 
constant in (2+l)-dimensional spacetimes p^ . 

In review of all these studies, now the following is clear: 

(a) In general the critical solution and the two dimen- 

sionless constants A and 7 are universal only with 
respect to the same matter field, and usually are 
matter-dependent. For example, for the collapse of 
the SU{2) Yang-Mills field, it was found H] that 
A K, 0.74 and 7 ~ 0.2, while in the case of massless 
scalar field, Choptuik found that A « 3.447 and 
7 « 0.37. 

(b) The critical solutions can have discrete self-simil arity 

(DSS) m or continuous self-similarity (CSS) 
or none of them, depending on the matter fields 
and regions of the initial data space. So far, in 
all the cases where the critical solution either has 
DSS or CSS, black holes form always starting with 
zero mass, and take the form of Eq.(2.7), the so- 
called Type II collapse, while in the cases in which 
the critical solution has neither DSS nor CSS, the 
formation always turns on with a mass gap, the so- 
called Type I collapse, corresponding, respectively, 
to the second- and first-order phase transitions in 
Statistical Mechanics |^ . 

(c) The universality of the critical solution and the ex- 

ponent 7 now are well understood in terms of per- 
turbations of critical solutions |l^] , while the one of 
A still remains somewhat of a mystery. The former 
is closely related to the fact that the perturbations 
of the critical solution has only one unstable mode. 
This property now is considered as the main crite- 
rion for a solution to be critical . 

To understand the last property better, let us con- 
sider the phase space, that is, consider GR as an infinite- 
dimensional dynamic system. If we make a (3 + 1) split 
of the spacetime, for example, following the Arnowitt, 
Deser, and Misner (ADM) decomposition, we will find 
that the dynamic quantities will be the induced spatial 
three metric, the extrinsic curvature, and the matter dis- 
tribution. Then, the phase space will consist of all the 
possible three metrics, extrinsic curvature, and configu- 
rations of the matter fields. For the case of massless 
scalar field, from the no-hair theorem of black holes 
we know that the only stable black hole solution of the 



Einstein-scalar field equations is the Schwarzschild black 
hole with a constant massless scalar field. Except for 
this black hole, another stable state is the Minkowskian 
spacetime. Of course, we also know that the collapse of 
a massless scalar field can form naked singularities, too, 
but so far we don't know if they are stable or not . At 
this point, we shall adopt the point of view of the cosmic 
censorship conjecture 0, and assume that they are not 
stable. Otherwise, there may exist two more critical solu- 
tions that separate, respectively, black holes from naked 
singularities, and flat spacetimes from naked singulari- 
ties. However, this doesn't affect our following discus- 
sions if we are restricted only to the boundary between 
black holes and flat spacetimes, and the analysis can be 
easily extended to other boundaries. 




Flat 

Space-Times 



FIG. 2. The phase space of the dynamic system of the Ein- 
stein-scalar field equations. The hypersurface S is the critical 
surface of codimension one, which separates the basin of black 
holes from the basin of flat spacetimes. A generic smooth 
family of initial data S\p\ always passes the two basins at the 
critical point p = p* on the critical hypersurface. All the 
initial data on the hypersurface will collapse to the critical 
solution A* that is a fixed point on the hypersurface when it 
has CSS and a fixed cycle when it has DSS. All details of ini- 
tial data are soon washed out during the collapsing process, 
and the collapse with initial data near the critical point will 
be very similar to the critical collapse. This similarity can be 
last almost to the fixed point A*, whereby the one unstable 
mode suddenly draws the collapse either to form black hole 
or a flat spacetime, depending on whether p > p* or p < p* . 



Then, we can see that the phase space can be divided into 
two attractive basins. Inside the dispersion basin, the flat 
spacetimes with different constant values of the scalar 
field, are attractive fixed points, while inside the basin of 
black holes, the Schwarzschild black holes with different 
masses are the attracting fixed points. The boundary 
between the two attractive basins is called the critical 
surface, and the critical solution always lies on it. Since 
it has only one unstable mode, this surface must be a 
hypersurface of codimension one, that is, one dimension 
less than the original infinite-dimensional phase space. 
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By definition, a phase space trajectory never leaves this 
hypersurface, if it is initially on it, but approaches to the 
critical solution, which is a fixed point on this hypersur- 
face if the critical solution is continuous self-similar, or a 
fixed cycle if the critical solution is discrete self-similar 
Ip^ . Within the complete phase space, the critical so- 
lution is an attractor of codimension one, i.e., it has an 
infinite number of decaying perturbation modes tangen- 
tial to the critical hypersurface and a single growing mode 
perpendicular to the hypersurface. Any trajectory begin- 
ning near the critical hypersurface, but not necessarily 
near the critical point (or cycle in the DSS case), first 
moves parallel to the hypersurface and goes down almost 
to the critical point (or cycle), then is suddenly drawn 
away by the single unstable mode in the perpendicular 
direction, and finally ends up at one of the fixed points, 
by either forming a black hole or a flat spacetime. Dur- 
ing the dynamic process, all details of the initial data 
are quickly washed away, except for the distance from 
the black hole threshold. Therefore, for the both super- 
critical {p > p*) and sub-critical {p < p*) collapse, there 
exists a domain, p* — Ap < p < p* + Ap, in the phase 
space, in which the collapse is very similar to the critical 
one during certain period of times [See Fig. 2]. 



III. CRITICAL COLLAPSE OF MASSLESS 
SCALAR FIELD: AN ANALYTIC TOY MODEL 



but slig htly different from the ones used in [Q. From 
Eq. ( |3.2|) it can be easily shown that the local mass func- 
tion is given by 



"^(w,w) = 2 



(1 - 4&2)uw 



(3.4) 



which is zero on the hypersurface v — and negative for 
u, V < 0. Thus, to have a physically reasonable space- 
time we need to restrict the above solutions valid only in 
the region u < 0, v > 0. Since the mass is zero on the hy- 
persurface V = 0, we may join the above solutions across 
the hypersurface v — with a Minkowskian spacetime. 
As shown in , this is possible if the metric in the re- 
gion V <0 takes the form of Eq. ( p.l| ) but with the metric 
coefficients and the massless scalar field being given by 
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r(u, v) — a{v) — —u — a(0), 
G{u,r) = 2a'{v), = 0, {v < 0), 



(3.5) 



where a{v) is an arbitrary function subject to a'{v) > 
and a'(0) = 1/2, and a prime denotes the ordinary dif- 
ferentiation with respect to the indicated argument. For 
such a matching, it can be shown that the hypersurface 
?; = is free of any kind of matter and represents a 
boundary surface Q]. The region w < is Minkowskian 
[See Fig.3]. 



In this section, we shall present a class of analytic solu- 
tions of the Einstein-scalar field equations, which repre- 
sents gravitational collapse of a wave packet consisting of 
massless scalar field |3^ . This class of solutions was first 
discovered by Roberts |4^] and later studied by several 
authors in the context of critical collapse [^ . As we shall 
show below, these solutions possess most of the features 
of critical phenomena, although they don't exactly repre- 
sent critical collapse, because the solution that separates 
the formation of black holes from that of flat spacetimes 
has more than one unstable mode [^. It is exactly in 
this sense, we refer these solutions as representing a "toy" 
model of critical collapse. 

The Roberts solutions are given by [M 



-G{u, v)dudv + r'^{u, u)c?^f2, 



(3.1) 



where u and v represent two null coordinates, in terms 
of which the metric coefficients and the corresponding 
massless scalar field 4> are given, respectively, by 



1 



Giu,v) = 1, 



(u-v)- (l- 462)i/2w 



{U - V) + {1 - ib2y^^v 



(3.2) 
(3.3) 



where 62 is an arbitrary constant. Note that the no- 
tations used here closely follow the ones used in |39] 




FIG. 3. The spacetime in the (u, u)-plane. The region 
u < is Minkowskian, while the region it < 0, v > rep- 
resents gravitational collapse of the massless scalar field, (a) 
When 62 < 0, the scalar wave collapses into a spacetime sin- 
gularity at u — —[(1 — 462)^''^ — l]w < 0, which is preceded by 
an apparent horizon located at it = — 4|&2|w < 0. (b) When 
b2 — 0, the spacetime singularity coincides with the apparent 
horizon on u = which is null, (c) When < 62 < 1/4, the 
massless scalar field first collapses and then disperses into in- 
finity, and finally leaves a Minkowskian spacetime behind in 
the region u, v > 0. 



On the other hand, from Eqs.(3.2) and (3.S) it can be 
also shown that the spacetime in the region u < 0,v > 
represents a collapsing massless scalar wave. When 
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b2 < 0, the scalar wave collapses into a spacetime sin- 
gularity on the hypersurface u = —[(1 — 4&2)^^^ — 
which is preceded by an apparent horizon at m = Ab2V. 
Thus, the corresponding solutions represent the forma- 
tion of black holes. When 62 = 0, the singularity coin- 
cides with the apparent horizon on the null hypersurface 
u — 0. When < 62 < 1/4, it can be shown that the 
massless scalar field first collapses and then disperses into 
infinity, without forming black holes, but instead, leaves 
a Minkowskian spacetime behind, which now is repre- 
sented by the re gion , u,v > 0, in which the metric takes 
the form of Eq.(|3.l[), but with 



G{u,v) =4bl^^b'{u), 



r = 62^%-5(u) + 6(0). 



1 + (1 - 4&2)i/' 



1- (1-462)1/2. 
(0 < 62 < 1/4, u,v> 0), 



(3.6) 



where b(u) is an arbitrary function, subject to b'(u) > 

1 /2 

0,6'(0) = 1/(4^2 ). One can show that the hypersur- 
face M = 0, w > is also free of any kind of matter and 
represents a boundary surface. 

I n th e case 62 < 0, where black holes are formed, 
Eq.(3.4) shows that on the apparent horizon u — ~A\b2\v 
the mass becomes unbounded as f — + -f-00. In order to 
have black holes with finite mass, we shall follow ||39| ] 
first to cut the spacetime along the hypersurface v — vq 
and then join the region < v < vq with an asymptoti- 
cally flat region. To model the out-going radiation of the 
massless scalar field, we shall choose the region v > vq a,s 
described by the Vaidya solution , 



(v > Vq), 



(3.7) 



where U is the Eddington retarded tim e, w hich is in gen- 
eral the function of u appearing in Eq.(3.1), and m{U) is 
the local mass of the out- going Vaidya dust. The corre- 
sponding energy-momentum tensor is given by 



where p is the integration constant, and 

, (1- 462)1/2 
P = z va- 



(3.11) 



For the details, we refer readers to Since for the 

above matching, the hypersurface v = vq is free of mat- 
ter, the function M{r) represents the total mass of the 
collapsing wave packet fill ed in the region < v < vq- At 
the past null infinity, Eq.(3.1C) shows that 



M{r 



-00) =p, 



(3.12) 



that is, the parameter p in the present case represents 
the total initial mass of the massless scalar wave packet 
with which it starts to c ollap se. 

As r ^ 0+, from Eq.( 3.10| ) we can see that Af (r) be- 
haves as 



M(r 



+00, 
0, 

-00, 



p> p 
p = p 
p<p 



(3.13) 



M(r) 



f = r/2 




FIG. 4. The mass M(r) of collapsing spherical shell. The 
line (a) corresponds to the case where p > p* , in which a black 
hole is formed, and its mass is given by Mbh = M(rAH)- The 
line (b) corresponds to the case where p = p* , while the line 
(c) corresponds to the case where p < p* . 



. 2 drniU) jj jj 



(3.8) 



The hypersurface v — vq in the coordinates {a;^} — 
{[/, r, 9, ip} is given by 



dUir) 
dr 



2r 



2m{U) 



, [v > Vq)- 



(3.9) 



Then, it can be shown that the junction conditions on 
the hypersurface v — vq require 



M{r)^m{U)l__ 



1 



p(4p*2 + r2)i/2 -2p* 



Ap* 



(3.10) 



On the other hand, it is well-known that the apparent 
horizon at r = 2M(r) of the out-going Vaidya solution 
always coincides with its future event horizon. Thus, by 
comparing the mass M{r) with r/2 we can tell whether 
the collapse forms a black hole or not. 



M(r) - - = 



4p* 



1/2 



Ar^[p+ {Ap* +r^)] 



(3.14) 



Clearly, only when p > p* , the scalar field and the null 
shell will collapse inside the event horizon at 



TAH = 2(p2 _ p* 



Nl/2 



(3.15) 
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FIG. 5. The corresponding Penrose diagrams, (a) The case 
where p > p* , in which a black hole is formed, and its mass 
is given by Mbh ~ M{rAH)- (b) It corresponds to the case 
where p = p* , in which the spacetime singularity becomes null 
and coincides with the apparent horizon, (c) It corresponds 
to the case where p < p* , in which no spacetime singularity is 
formed, instead, when the wave packet of the massless scalar 
field collapses to the hypersurface u = 0, all of its mass has 
been radiated away and nothing is left to collapse, so the 
spacetime in the region u, v > becomes Minkowskian. 



When p — p*, M{r) — r/2 is possible only at the origin, 
r = 0, where a zero-mass singularity is formed. Thus, 
the solution with p = p* represents the "critical" solution 
that separates the supercritical solutions {p > p*) from 
the subcritical ones {p < p*). In the subcritical case, 
M{r) is always less than r/2, and the collapse never forms 
a black hole [cf. Fig. 4]. 

In the subcritical case, the region u,v > should 
be replaced by the Minkowskian solution ( ^.6[ ). As 
shown above, the matching across the hypersurface u = 
0, < u < Wo is smooth, i.e., no matter appears on 
it. To show that it is also the case on the hypersu rface 
M = 0, w > Wo, which separates the Vaidya solution (3.7) 
from the Minkowskian one (^^), we first make the co- 
ordinate transformation U — U(u), and then write the 
metric (3.7) in terms of u. Using the results obtained in 
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p6[ , one can show that to have a smooth matching we 
have to impose the condition 



U'iO) = 



1 



(462)1/2' 

M(r) |„=o = . 



(3.16) 
(3.17) 



Clearly, by properly choosi ng th e function-dependence of 
U(u), the first condition ( 3.16 ) can be always satisfied. 
On the other hand, from Eqs.( |3.2| ) and (3.10) one can 
show that the last condition is also satisfied identically. 
Therefore, the matching of the Vaidya solution to the 
Minkowskian one across the hypersurface w = 0, v > vq 
is always possible for p < p* . The corresponding Penrose 
diagram for each of the three cases are shown in Fig. 5. 

On the apparent horizon r = tah, the total mass of 
the scalar wave packet is given by 



MBH^lrAH = K{p~p*)^/^, 



(3.18) 



where K = (p + p*)i/^. The above expression shows 
that the black hole mass takes a power-law form with 
its exponent 7 = 0.5, which is different from the DSS 
case, where Choptuik found 7 « 0.37. As we mentioned 
previously, the above solutions don't really represent crit- 
ical collapse, because the "critical" solution given above 
has more than one unstable modes |^^. Therefore, the 
different value of 7 obtained here does not means any 
contridiction to Choptuik's numerical results. 



IV. CURRENT INVESTIGATIONS 

As Critical Phenomena in gravitational collapse is a 
rather new area in GR, there are many open problems. 
In the following we shall mention some of them. 



A. The Effects of Angular Momentum in Critical 
Collapse 

As we know, the angular momentum plays a very sig- 
nificant role in black hole physics, and all the realistic 
bodies, such as, neutron stars, have non-zero angular 
momentum. Thus, it is very important to study the ef- 
fects of angular momentum on critical collapse. So far, 
all the studies of critical phenomena in gravitational col- 
lapse have been restricted to spherical case, except for the 
works of Abrahams and Evans ||l^ and Alcubierre et al. 

In jl^ the authors studied the collapse of axisym- 
metric purely gravitational waves, and found the type 11 
critical collapse. However, in this study the total angular 
momentum is still zero. In the collapse of pure Brill 
type gravitational waves in 3D Numerical Relativity was 
studied and the critical amplitude for black hole forma- 
tion was determined. However, due to the complexity of 



the problem, no sufficient evidence for critical collapse 
was observed. 

In addition to the above, Gundlach and his co-workers 
psi studied the problem using non-sphcrical perturba- 
tions, and in particular found that all the modes of non- 
spherical perturbations of the massless scalar field stud- 
ied initially by Chpotuik are stable, and, as a result, 
Chpotuik's critical solution found in the spherical case 
may remain critical even in non-spherical case. He also 
found that small angular momentum also takes a scaling 
form near the critical point but with a different expo- 
nent. Besides, Rein, Rendall and Schaeffer studied 
the spherical collapse of coUisionless matter that consists 
of counter-rotating particles, and found that only Type 
1 critical collapse. This result was further confirmed by 
Olabarrieta quite recently . 

Moving from spherically symmetric case to axisymmet- 
ric one, the problem becomes much mathematically in- 
volved, and very sophisticated (numerical) methods are 
needed. Choptuik, Hirschmann and Liebling, among oth- 
ers, have been working on this problem recently W9|, and 
are expected to report their results soon. 



B. The Quantum Effects on Critical Collapse 

Critical phenomena are actually phenomena in the 
strong gravitational field regime, and Quantum effects 
should be very important for the formation of black holes 
with very small mass. Chiba and Siino studied this 
problem and showed that the Quantum effects may de- 
stroy the type 11 critical phenomena, while Ayal and Pi- 
ran |5l| showed that they don't, but rather shift the crit- 
ical value p*. However, since in both of the two cases the 
Quantum effective energy-momentum tensor (EMT) was 
taken from two-dimensional toy model, the consistence 
of such an EMT with the four-dimensional gravitational 
collapse is still an open question. Recently, Brady and 
Ottewill [ID calculated the effective EMT of a confor- 
mally coupled scalar field on the fixed background of the 
critical solutions of the perfect fluid with the equation 
of state p — kp in four-dimensional spherical spacetimes, 
and found that when k < 0.53, the Quantum effects de- 
stroy the type 11 critical phenomena, while when k > 0.53 
their calculations break down, and a definitive conclusion 
is still absent. 



C. The Application of Renormalization Group 
Theory to Critical Collapse 

The Renormalization Group Theory has achieved great 
success in the studies of critical phenomena in Statistical 
Mechanics , and several authors, including Argyres 
and Koike, Hara and Adachi ||l8|, have pointed out 
that the time evolution near the critical solution in grav- 
itational collapse may also be considered as a renormal- 
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ization group flow on the phase space of initial data. As 
a matter of fact, the analysis of the phase space given in 
the Introduction exactly followed this idea. However, this 
analysis is valid only in self-similar spacetimes. In order 
to obtain a full rcnormalization group, one needs to gen- 
eralize them to arbitrary spacetimes, which is turned out 
not trivial, as in GR the choice of coordinate systems is 
completely arbitrary. Garfinkle and Gundlach |]5^ have 
taken some initial steps to this direction, but a success- 
ful application of the Renormalization Group Theory to 
critical collapse still remains as an open question. 

Besides the above mentioned problems many others 
are also under the current investigations, such as, finding 
more matter fields that exhibit critical collapse, includ- 
ing universal classes; applying the analysis of perturba- 
tions of black holes to critical solutions; understanding its 
physical origin of the constant A |55|; finding some possi- 
ble astrophysical observations of critical phenomena, and 
so on. In particular, it was known for a long time that 
the collapse of neutron stars exhibits the type I critical 
phenomena |5^]. An important question is that: Does 
this type I critical collapse have any observational conse- 
quence? 

For further references of critical phenomena in gravita- 
tional collapse, we would like to refer the readers to the 
review articles kHISI- 
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